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Abstract 

We investigate optimal control strategies for state to state transitions in a model of a quantum 
dot molecule containing two active strongly interacting electrons. The Schrodinger equation is 
solved nonperturbatively in conjunction with several quantum control strategies. This results in 
optimized electric pulses in the THz regime which can populate combinations of states with very 
short transition times. The speedup compared to intuitively constructed pulses is an order of mag- 
nitude. We furthermore make use of optimized pulse control in the simulation of an experimental 
preparation of the molecular quantum dot system. It is shown that exclusive population of certain 
excited states leads to a complete suppression of spin dephasing, as was indicated in Nepstad et 
al. [Phys. Rev. B 77, 125315 (2008)]. 



PACS numbers: 73.21.La, 78.67.-n, 85.35.Be, 78.20.Bh 



I. INTRODUCTION 



Application of quantum control theory to optimize transitions in strongly interacting 
quantum systems is a well established technology in simple two-level systems pQ. In more 
complex or open systems involving many important channels it is in general much more 
complicated to improve the probability of desired reactions and transitions. This is true 
in systems as diverse as dipole blockade dynamics in cold Rydberg gases [2J and electron 
dynamics in semiconductor two-electron quantum dot systems [3j. In the latter, which is our 
case, the electron-electron interaction is comparable to other interactions in the system and 
cannot be neglected. The ability to achieve fast and optimized transitions in such systems, 
and a variety of others, is important for improving present day technology in quantum 
information, metrology and in quantum chemistry. 

In few-electron quantum dots it is well recognized that interactions with the substrate 
will induce decoherence, which limits the ability to utilize unique quantum properties such 
as entanglement. Examples of such interactions are hyperfine and spin-orbit interactions 
between the quantum dot electrons and surrounding atoms, and interactions with phonons 
in the substrate lattice. As strategies to reduce decoherence, one can either carry out 
experiments in systems and at temperatures which minimize unwanted interactions, or try 
to develop methods to perform the required transitions much faster than the characteristic 
timescale of the decoherence. We have previously demonstrated that intuitively selected 
microwave pulses can populate both single states and more complex states of the lowest 
excitation bands, and we were able to further decrease the transition time in the first case 
by optimal pulse control. [3] 

In the present work we optimize time- dependent transitions to more complex target 
states and compare various strategies of optimization including frequency-selective control 
algorithms. [3] We show that more advanced control strategies lead to a factor of 7 faster 
transition times than previously reported using intuitively constructed pulses. In the second 
part we address the application of quantum control inside regions of anticrossings. This 
is related to a recent experiment by Petta et ai, [5] which measures spin dephasing of 
the system through hyperfine interactions with the surrounding nuclear spin bath. The 
experiment was simulated in Nepstad et ai, [6J and very good agreement between theory 
and experiment was achieved. In the same work, we further demonstrated how populating 
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higher excited states could be used as a method to inhibit decoherence. In this paper we 
apply the technique of optimal control theory to exclusively populate such states during 
initial setup of the experiment. The following section describes the theory in detail. In 
section III we present the results, followed by concluding remarks. 



II. THEORY AND METHOD 



In this section we review and detail the numerical methods used to study dynamics of 
a two-dimensional, two-electron double dot exposed to electric and magnetic fields. This 
includes DC and pulsed electric fields, strong external magnetic fields and weak, locally 
varying magnetic fields representing the hyperfine interaction . [31 16|, [7j 



A. Model 



The two-dimensional single-particle effective mass Hamiltonian of our system reads 
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Combined with the electron-electron interaction term, the total field-free Hamiltonian be- 
comes 

Ho = /lo(ri) + hoir^) + . (2) 

In addition, we include external magnetic and electric fields, 
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+ -feBextSz - e^{t)x, (3) 

and define the total field Hamiltonian as Hext = hext{T^i) + ^ext(i"2)- Here ri^2 are single- 
particle coordinates, ^ is an electric time dependent field applied along the inter-dot axis 
and Bext is an external magnetic field perpendicular to the dot. The material parameters 
may take on different values to reflect various physical systems. In this paper we will use 
values compatible with GaAs quantum dots, where m* = 0.067me (effective mass), = 12.4 
(relative permittivity), 7e = 9* 2^ (gyromagnetic ratio) and g* = —0.44 (effective g-factor). 
The electron mass is denoted mg. The confinement strength is set to hu = 1 meV and the 
interdot separation to d = 130 nm, which are realistic experimental values. [HI |9] 



We obtain eigenstates of the field-free two-electron Hamiltonian Hq by Arnoldi 
iterations |lUj using a basis of symmetrized products of one-electron harmonic oscillator func- 
tions, 0„ ^ = E">r Cijl^j)^ where 

M = < ^ (4) 
[0i(ri)0j(r2) « = j, 

with i, j = {nx,ny} representing the quantum numbers of the Hermite polynomials in x and 
y respectively. The symmetric and antisymmetric basis functions correspond to singlet and 
triplet states respectively. We obtain converged results in all cases using Uy^max = 4 and 



B. Dynamics in the eigenstate basis 

In cases where the hyperfine interaction between the two active electrons and the semi- 
conductor nuclei surrounding the quantum dot can be neglected, the total spin is a conserved 
quantity. We then need only consider the subspace of symmetric basis functions, correspond- 
ing to singlet states, choosing the + sign in Eq. Q. The dynamics is governed by the time 
evolution of the expansion coefficients, 

^c^j{t) = Y^Ci,j,{t){iljl\H\ij). (5) 
iiji 

This system of equations is then integrated using an adaptive form of Adam's method. [TT] 
In the singlet subspace using a basis of ~ 4000 states, the calculations are reasonably fast. 
A considerable speedup may be obtained by switching to the eigenstate basis. In this case 
propagation times of nanosecond duration is performed in less than a minute (on a dual core 
AMD Turion 64 bits processor). We find converged results using a basis of 50 eigenstates. 
The coefficients in Eq. (|5]) become the coefficients of the eigenstates 

idk (t) = J2di{t){l\ Hext \k)+Ekdk{t), (6) 

where now |/) indicates eigenstate / with corresponding energy Ei. The corresponding field 
matrix elements are calculated using analytic expressions obtained in the harmonic oscillator 
basis. The required matrix elements are given in detail in Popsueva et al. [7] 



4 



C. Hyperfine interactions 



A particular source of decoherence in double quantum dot molecules is the hyperfine 
interaction with the surrounding substrate nuclei, which has a characteristic timescale of a 
few nanoseconds. [I2] To study this interaction from first principles spin couplings for ~ 10^ 
nuclear spins surrounding the electrons must be included. The magnitude of the interaction 
is consistent with a random magnetic field of a few mT. For the timescale of the experiment 
( 50ns) this is taken to be constant and its spatial dependence can to a good approximation 
be represented by a step function, [B] 

{{BxGx + ByGy + BzGz) , for X > 
^ ' ' ^ (7) 

, otherwise. 

The hyperfine interaction term is then given by 

H^ = -1,Y,S,- Ba,, (8) 

i=l,2 

where Sj is the spin operator of electron i. In this semiclassical picture, we must con- 
sider an ensemble of quantum dot systems, each with a different random nuclear mag- 
netic field, and average over the ensemble to obtain physical quantities. To obtain 
the ensemble, we use a normal distribution of magnetic fields about zero, P{Bn) = 
l/(27ri?^^j.)5 exp (— -Btv • /2B'^^^^). [13] i?nuc can be determined by experiments and is of 
the order of 1 mT. [9j The interaction term induces couplings between the singlet and triplet 
states and between the different triplet states, necessitating the inclusion of both subspaces 
in the calculations. Details of the matrix elements involved can be found in Appendix A. We 
remark that other types of interactions with external degrees of freedom, such as interactions 
with electron spins or phonons, can be introduced formally in the same way. 

D. Dynamics in the adiabatic basis 

Experimental studies of spin dephasing in quantum dots require preparation of two elec- 
trons in the singlet ground state. In a recently reported experiment by Petta et ai, [5] this 
was achieved by applying a large external electric field to the double dot, deforming the con- 
fining potential until at large field strength it became essentially a single dot, as illustrated 
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FIG. 1: Illustration of a spin dephasing experiment in a double quantum dot molecule, as we model 
it in our simulations. The blue (dark) line shows the slowly switched electric field that guides the 
two electrons from a single- to a double-well configuration. During the switch control pulses may 
be applied in order to guide the electrons into excited states. 

in Fig. [T] The electrons were then allowed to tunnel into the trapping region, forming a sin- 
glet state. Reversing the electric field slowly guided the electrons into the ground state of a 
delocalized double-well configuration, where dephasing occurs. Final readout was performed 
by once more tuning the electric field to a single dot configuration. 

To simulate such an experiment, the one-center basis approach as described above is 
unsuitable, as a very large number of basis states is required to accurately represent the 
wavefunction when the electric field is large. Including the triplet states adds an additional 
factor of four to the basis size, making the calculations prohibitively time consuming. Even 
switching to the diabatic b described above, yields lengthy calculations and conver- 

gence in terms of basis size is arduously obtained. However, we observe that the energy 
spectrum as a function of electric field strength displays well-separated states with clear 
anticrossings. These considerations lead us to consider instead an adiabatic basis approach, 
where the wavefunction is expanded in eigenstates depending parametrically on the electric 
field, 

^(n, r2, t) = J2 Ck{t)ek{r,, ra; ® l^) , (9) 

k 

where IS") refers to either a symmetric (triplet) or antisymmetric (singlet) spin function. 
Note that the electric field is time-dependent {C, = ^(t)), but we have dropped the explicit 
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reference to t in order to simplify notation. The basis states 6^ are determined from the 
eigenvalue equation 

(Ho - eiX) ^(n, rs; = e(e)^(ri, v,- 0- (10) 



Inserting Eq. ^ into the TDSE and using Eq. 10, we find the governing equation for the 
coefficients. 

Ck{t) = -ei J2 ^^^^^^c,(t) + ze,{Ock{t), (11) 

Written more compactly on vector form, this reads 

c(t)= (-eeK(e)+2£(0)c(t). (12) 

The antihermitian matrix K(,^) is computed for a set of electric field values {^m} using the 
numerically obtained eigenstates and eigenvalues together with analytic matrix elements of 
the symmetrized harmonic oscillator functions defined in Eq. |4| 

where index m refers to the electric field points. Since the explicit time dependence in 



Eq. (12) is only found in the scalar function ^(t), the matrix elements need only be computed 
once, speeding up the time integration. Only ^(t) must be computed during integration, but 
this is inexpensive. As the numerically computed basis set is not continuous in ^, we use a 
simple interpolation between the ^ grid points where required by the integrator. 



E. Optimal control 

In this section we describe the iterative Krotov method [H] for optimizing optical tran- 
sitions in quantum systems. [I5] Following we will describe and discuss modifications to 
this scheme. In general the method aims to maximize the expectation value of a positive 
semi-definite operator by means of an external field while minimizing the field energy. The 
time evolution of the system in which we want to optimize transitions is described by the 
time-dependent Schrodinger equation 

z^^^{r,,r2,t) = [Ho-ee{t)X]^{r,,r,,t), (14) 
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where e(t) is an electric field and X = xi + X2 as before. We have chosen to use e(t) for 
the electric field whenever we refer to (optimized) pulses. Hq is the field free Hamilto- 
nian in Eq. Our goal is to apply control theory to find optimal pulses for population 
transfer from an initial state $j = \l'(t = 0) to a target state $4. The states are prese- 
lected and the pulse duration is fixed to t = T. The optimization is done by maximizing 
the expectation value of a projection operator, that is, maximizing the functional 

J^[^] = (^f(T)|$t)(<l>t|^(T)) = |(<I><|^(T))|2. The requirement that the field intensity should 
be as small as possible is achieved by minimizing a second functional, J2 [e] = dt A(t) [e^ (t)] , 
where the predefined function A(t) acts as a penalty factor, which can be used to impose 
an envelope on the electric field. We will use X{t) = A unless otherwise stated. In each 
iteration, the updated control field is found as a solution to 

VJ„[e]=0, Ja = Ji-J2. (15) 

We proceed to sketch a simple implementation of the Krotov iteration algorithm: The time 
interval, [0,T], is divided into fixed-length intervals ti on which e{t) is taken to be constant. 



e{ti) = et^, ti G [0,T]. The first step is to integrate the initial value problem of Eq. (14). 
For the first iteration, 1 = 0, use some initial guess for the control, e°(tj). The choice of 
initial control is by no means immaterial, as we will see later. After propagating forward, 
calculate the yield, |($t|\l/-'^(T))p, where \1/'^(T) is the final state. If the desired yield has 
been reached, the iterations are terminated. If not, solve the terminal value problem, 

X = -^ [Ho - fie{t)] x, with x(T) = \<^t) {<l>t\^\T)) , (16) 

and obtain x\^)- The updated control components ef^^ are obtained while integrating the 



TDSE, Eq. (14), again (step-wise) with \I'^+^(0) = For the first time interval choose 
e[+^ = -Im(x^(0)|/i|*^+^(0))/A, and with this e[+^ integrate to find ^^+^(ti). Repeat the 
process for the next time interval, using 



The entire procedure is repeated until maximum number of iterations or desired yield is 
reached. Using the eigenstates as basis makes it possible to perform several hundred itera- 
tions in a few hours. 



The expression for ej., Eq. (17), is a zeroth order approximation to the update equation 



given to full order in Degani et al. [1] Although the simple iteration method described above is 
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not guaranteed to converge monotonically, it works quite well for the system at hand. In fact 
this feature might even be desirable, as it acts as a 'shake-up' of the numerical calculations: 
the iteration scheme finds only local maxima and thus adding small perturbations to the 
solution through having a non-monotonous convergence might lead to even better optimal 
controls. Indeed, this effect was observed when performing an additional update of the 



control during the backward integration, cf. Eq. (16). In this case the convergence is 
smoother but often stagnates around a lower maximum yield. However, we would like to 
point out that the technique of using backward updates has proven to be quite effective in 
simpler systems, acquiring extremely high yields. IH [H] Tests using a basis of only a few 
states confirmed this also in our system. 

The method as presented above has another restriction in that it does not discriminate 
between possible controls, except favoring those of low intensity. This often leads to quite 
complicated controls that are difficult to produce in an experimental setup. There have been 
some attempts to address this problem. [T71 [TH] In Werschnik and Gross, a desired structure 
is enforced by projecting the control onto a preferred subspace in every iteration. Instead of 
this brute force strategy, a modified functional, J{,, can be introduced |1] which favors low 
energy controls with a desired structure. This is achieved by choosing a set of 'good' controls 
spanning a desired subspace of the full control space. The 'bad' control subspace is then 
defined to be the orthogonal complement of the 'good' subspace. The weighted terms of the 
projection of the control onto the 'good' and 'bad' subspaces are added to the functional J^. 
The new functional to maximize is 

Ji[e] = (*(r)|$t)($i|^(T)) 

- dtU^{\I + Xi^good + M^bad)U, (18) 

where the A-term is J2, and ^good/bad are projection operators onto the 'good' and 'bad' 
control subspaces respectively. Here we have used a set of frequencies corresponding to 
transitions between the 10 lowest bound states as our space of 'good' controls. More specif- 
ically the space of 'good' controls is defined as the span of /jj, z = 0, . . . , 9, j < i, 

fi.itk) = sm\7itk/T) cos {{Ei - Ej)tk) , (19) 

t € [tk,tk+i]- [IH] Optimizing Jb instead of Ja guides the control algorithm in the direction 
of controls with desired frequencies. 



III. RESULTS 



We here present results of calculations based on the control schemes outlined in the 
previous section, with respect to optimizing simple and combined state to state transitions 



in the double dot system. Section |III A| deals with transitions in the singlet subspace at zero 
electric field, while section III B focus on transitions during electric gate switching, indicated 
in Fig. [2] at the point of anticrossing (dashed black circle). 

The eigenvalue spectrum in Fig. [2] is shown as a function of electric field strength (left). 
As an initial strong negative electric field is decreased the state energies are seen to increase 
linearly and a number of anticrossing regions appear. The physics of the anticrossings 
normally involves strong state mixing. For example, the ground state in the circled area 
changes from a state containing essentially both electrons in one dot, to a covalent two-center 
state as the anticrossing is traversed. The molecular states at zero electric field (right panel) 
were classified and labeled in Popsueva et al. [7j 

In the present calculations we employ electric fields linearly polarized in the x-direction. 
This couples states which have different a;-parity and equal y-parity. Fig. [2] shows only states 
with the same y-parity as the ground state. The states |1), |2) correspond to single exciton 
states while the states |5), |6) are ionic states. As the dot separation tends to infinity, these 
become degenerate and consist of two-electron single dot ground states with both electrons 
in the same dot. 
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FIG. 2: The two-electron double dot spectrum as a function of electric field strength (left), and 
details of the spectrum at zero electric field (right). Solid lines (blue) are singlet states, dashed 
lines (red) are triplet states. The arrows indicate transitions referred to in the text. Only states 
of even y-parity are shown. 
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A. Spin conserving dynamics in the singlet subspace 

First, neglecting any spin interactions, we restrict our attention to dynamics in the sub- 
space of singlet states. From the ground state, transitions to the states |2) and |6) are 
dipole-allowed, while state |5) can be reached via |2). We will study each of these tran- 
sitions, finding that optimization procedures can produce very short pulses which achieve 
almost unit probability transfer. 

1. Single state transitions 

Investigating the |0) — |2) transition, we find that an "intuitive" sine-squared envelope 
pulse tuned to the resonance frequency will transfer 98.7% of the population in 237 ps. The 
population of |2) during the pulse is shown in Fig. |3| labeled I (black curve). In a previous 
attempt at optimizing this transition, [3J we found that 96.5% transfer could be achieved 
in 111 ps, using an energy penalty functional and amplitude cutoff (II - gray curve). With 
the present approach, the same functional (J^) provides better results, transferring 98.6% 
of the population in only 67 ps (III - blue curve). Replacing the energy functional Ja with 
the structure functional Ji, gives a slightly better final population of 99.3% (IV - red curve). 

With transition time decreased to 67 ps, the population transfer proceeds in an irregular 
manner for the energy-penalty optimized pulse. III. During the pulse, as much at 70% of the 
population is transfered to highly excited states (> 10). Direct transitions to these states 
from the ground state may be discouraged by using the structure penalty which will favor 
the corresponding frequencies and disfavor all others. Switching to the functional Jb, the 
resulting pulse causes population of higher excited states to reduce to 20% (IV) . The high 
numbers are mainly due to the 13th excited state, which has a strong coupling to the 2nd 
excited state. The resonance frequency of this transition is close to that of the |0) — |2) 
transition, which is included in the space of good controls. Disregarding the population of 
1 13), the population of the remaining higher excited states is 31% and 2.7% for the Ja and 
Jh optimizations respectively. The plateau structure in the population of |2) during the two 
short pulses (III and IV) is due to transient population of |13). 
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FIG. 3: Figure showing (optimized) transitions from |0) to |2). Black line (I): "intuitive pulse", 
a{t) cos{iot) , where lu = {E2 — E^j/h = 1.5 THz and a{t) is a sin'^-envelope. Gray line (II): 
optimized pulse using the functional Ja, duration is 111 ps. Red line (III): optimized pulse using 
the functional Ja, duration is 67 ps. Blue line (IV): optimized pulse using the functional Jj,, 
duration is 67 ps. The optimization was done using At = 0.5 

2. Charge localization 

Previously we demonstrated how charge localization in one dot can be achieved in less 
than a nanosecond by applying weak, resonant pulses on the system. [3] The charge localized 
state (CLS) is a combination of two states in the third energy band of the spectrum exhibiting 
ionic structure, in analogy to ionic states in diatomic molecules. At large interdot separation 

the two states resemble the asymptotic states 

\g{riL, r2L)) ± \g{riR, ra/?)) (20) 

where \g) refers to the shifted ground state of a single two-electron dot. Creating an equal 
linear combination of these states will cause the two electrons to oscillate between localization 
in the left and the right dot with a period of 180 ps, inducing a weak current over the 
dot. This is illustrated in Fig. |4], where the upper panel shows the expectation value of 
X = xi + X2 as a function of time, its value oscillating between the two minima of the 
double-dot potential. Also shown is the integrated one-electron density of the CLS at certain 
times during field-free time evolution, p(x) = / dyi d'^r2 |\l/(ri, ra)^. 

In the intuitive scheme the transition to the CLS is achieved via an intermediate transition 
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FIG. 4: Field-free time evolution of the charge localized state. Upper panel: Time evolution of 
the expectation value X = xi + X2- Lower panel: Single-electron density averaged over the y- 
coordinate at three different times during the time evolution, t = ps (left), t = 45 ps (center) and 
t = 90 ps (right). The red markers indicate the value of {X) at the three times. The confinement 
potential is indicated by the dashed lines. 



to the second excited state in the second energy band (labelled |2) in Fig. |2]). This is 
necessary because the lowest ionic state has positive x-parity and can not be reached from 
the ground state directly, due to selection rules. The three transitions involved are indicated 
by arrows in the rightmost part of Fig. |2} Fig. |5] (second panel from top) shows the eigenstate 
population as a function of time during the sequence of resonance pulses and during the 
optimized pulse (bottom panel). The respective pulses are shown above. The first two 
pulses in the uppermost panel use a sin^ envelope whereas the last pulse uses a sin^ ramp-on 
over 10 oscillations. The optimized pulse was obtained using the functional Ja with final 
time T = 117 ps and a maximum of 300 iterations. We used A(t) = 1/ sin^(7rt/T) to ensure 
that the pulse is zero at t = and t = T (note that A is a penalty factor, making the penalty 
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for a non-zero field at the endpoints infinite). As seen from the population during the 
optimized pulse, the strategy of using the second excited state as an intermediate transition 
is also being used here, although the transitions are somewhat more involved. The total 
transition time has been brought down from 852 ps to 117 ps using the optimized pulse, and 
the population of the target state has been improved from 97.2% to 99.8%. An important 
thing to note is that using a defined target state we are also able to selectively choose the 
configuration of the charge localized state which is determined by the relative phase between 
the ionic states. We have used a target state defining the two electrons in the left dot at the 
end of the pulse. The charge oscillations of the CLS has a period of ~ 180 ps, and so the 
relative phase evolution of the ionic states is important during the propagation. In this sense 
the optimal control scheme is stricter compared to the intuitive approach, where we did not 
control the final configuration of the electrons, only the population in each eigenstate. When 
considering simple state to state transitions, the phases naturally are not important. In 
Fig. |6] we have applied to optimize transition to the charge localized state. The results 
for T = 67 ps are compared with optimization using the Ja functional. Results for Jf, and 
J a are shown in blue and red curves respectively (all panels). The upper panel shows the 
convergence of the yield (projection of the final state onto the target state) as a function of 
iteration. Additional iterations did not produce higher yields. The red curve (upper panel) 
is fluctuating strongly and the final pulse is also somewhat irregular. In this case we applied 
a low pass filter to the final pulse to get rid of very high frequency components, caused by 
numerical noise. We checked that removal of these components did not affect the final yield 
and dynamics. In general, we experienced greater difficulties in achieving converging results 
using only the energy penalty, and the yield often converged to zero. The maximum yields 
for the two methods were 94.3% (Ja) and 83.3% {Jh)- While the structure penalty strongly 
limits the presence of unwanted frequencies in the optimized pulse, population of excited 
states beyond the 10 lowest still occurs. This is again related to the existence of resonant 
transitions to higher excited states matching the frequency of the desired transitions. The 
population of excited states (> 13) is for both methods ~ 20%. 

We have noticed in all our calculations that the pulse produced by the optimal control 
algorithm is sensitive to the choice of initial field. An example of this is shown in Fig. [7j Here, 
we start the iterations using a constant initial field (gray, horizontal curve), and consequently 
obtain a rather different optimal pulse compared with the one in Fig. [6] (blue curve), where a 
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sin^ enveloped pulse was used (gray curve). In this case we see that the optimized pulse has 
retained much of its initial DC component. The maximum population of the target state is 
98.6% after 419 iterations. In this case the population of highly excited states during the 
pulse is considerably less, with < 10% in the 24 highest states. Optimization using only the 
energy penalty in this case gave a very short, high frequency and high intensity pulse, with a 
resulting yield of only ~ 60%. These examples illustrate the limitations of using only energy 
penalty when the propagation time becomes short, and how adding structure penalty can 
consistently guide the control towards a wanted frequency space. 

We end this section with some comments on the issue of using y-polarized fields. The 
transition to the lowest ionic state (with positive a;-parity) could in principle have been 
achieved using a y-polarized field and the third excited state in the energy spectrum, which 
has positive x-parity and negative y-parity. Note that the transition to the upper ionic state 
(with negative x-parity) can only be reached via x-polarized fields. Using this scheme one 
could perform the two operations simultaneously with weak fields. 

There are however properties of the spectrum obstructing the use of y-polarization in this 
system. In the case of the CLS, the coupling between the third excited state and the lowest 
ionic state is virtually zero. Moreover, from the lower energy bands, there exist a multitude 
of strong couplings to states further up in the spectrum, precluding selective transitions 
to lower lying states. Numerical calculations confirm that selective state population is 
impractical using ?/-polarized fields. 

This feature of the double dot spectrum is related to the symmetry of the potential, 
particularly in the y-direction. As noted in an earlier work, [3] and as we will also see 
manifested later on, optical manipulation in this two-electron molecule system is actually 
restricted by the degree of symmetry in the potential, and control would be more easily 
achieved in slightly asymmetric dots. Similarly, we expect that ?/-polarized fields could be 
more useful in anharmonic systems. 



B. Optimized transitions and spin interactions 



In a previous paper, [Bj we studied the effects of spin-dephasing in the quantum dot 



system, modelling the experiment described in Sec. IIP and replicating experimental con- 
ditions as accurately as possible. We observed that when using ultrafast electric switching 
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(1 ps) through the anticrossing (black-dashed circle in Fig. |2]), large population transfer from 
the ground state to the second energy band resulted. The decohence is largely suppressed for 
those states, and when the system was switched back to the 'single dot' configuration, 95% 
of the initial singlet population was regained. The suppression is explained by the fact that 
at zero electric field, the singlet-triplet energy splitting is approximately 100 times greater 
for the second excited singlet state compared with the ground state. Some of the population 
vanished to higher excited states during passage through the anticrossings, causing the 5% 
loss. 

By applying optimal control schemes in combination with the adiabatic electric switch, 
the transition to excited states may be achieved with near 100% probability. An optimized 
pulse applied at the point of anticrossing will force a non-adiabatic transition and by target- 
ing the desired excited state explicitly we minimize loss to other states. Fig. |8] shows such 
a transition between the two lowest eigenstates using an optimized pulse. The optimized 
pulse was obtained using the structure functional Ji, together with \{t) = 1/ sm{nt/T), and 
has a duration of T = 67 ps. In this case the population of other states during the pulse 
is completely negligible and the final population of the second excited state is as high as 
99.9%. After the pulse, an adiabatic switch of duration 2 ns is applied, guiding the system 
to the delocalized double-well configuration, where the system is left to interact with the 
spin bath for 50 ns. Reversing the adiabatic switch and optimized pulse procedure, we find 
that 99.3% of the ground state population is regained. 

IV. SUMMARY 

In this work we have demonstrated to what extent quantum control strategies can be 
applied to obtain required transitions between electronic states of two-electron quantum 
dot molecules. Such transitions are non-trivial partly due to the strong electron-electron 
interaction, but also the large number of coupled states induced by the external fields. 
Nevertheless, the calculations have shown that single states and superposition of states may 
be reached with close to 100% probability. 

Using weak, pulsed electric fields in the THz regime, we have shown that transitions 
from the ground state to a preselected excited state may be obtained within 100 ps. When 
applying advanced control strategies, a speedup of more than 7 times the transition time 
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using straightforward intuitive pulses is gained. Such control strategies also have the ad- 
vantage of returning a final pulse consisting of experimentally relevant frequencies. In the 
case of interactions with slowly varying external fields, which have been applied in experi- 
ments, we have shown that complete transitions at anticrossings can be obtained. This is 
a realistic implementation of a fully diabatic time development, which in the Landau-Zener 
model requires infinitely fast transitions. We also showed that the hyperfine interaction in 
the excited states is unimportant at the considered time scales as opposed to in the ground 
state. Advanced engineering of tailored pulses as here described appear as a realistic route 
to accessing and manipulating electronic states in experiments. 



Appendix A 

Matrix elements for the hyperfine interaction 

In this Appendix we give details of the matrix elements for the hyperfine interaction in 
Eq. ^ using the symmetrized basis of Hermite functions. The spin-states are as usual 

' X^(l)xn2) = |T+) 

Triplet Ki(l);^i(2) = |r__) (21) 

^ ^(X^(1)XH2) + XK1)X^(2)) =|To) 

Singlet (xHl)xH2) - xKl)xH2)) = \S) (22) 
The singlet has corresponding symmetric spatial function of the form 

1*7(^,^2)) ^ <^ ^ (23) 

[ Kj) = 73 (0.(^1)0.(^2) + 0,(ri)0,(r2)) 

The triplet has corresponding antisymmetric spatial function of the form 

|^j(n, r^)) ^\kl) = ^ {Mri)Mr2) - Mri)Mr2)) (24) 

The 4'i{r)s are as before two dimensional Hermite functions with i = {nx,ny}. Recall also 
the representation of the effective nuclear field, Eq. ([T]), 

{B^e^ + ByGy + B^e^) , for x > 
Bat = \ (25) 

, otherwise. 
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Matrix elements for the \S) < — > |T) coupling 



(^,(n,r2);^|^^|^j(ri,r2);T) (26) 



The first case reads 



{^^■,S\J2S^■BM\kl■,T) = 

i=l,2 

J2 {Su{^\B.,^\k)-5,k{^\B,,^\l)} (27) 

^ Xi=x,y,z 

X {S\ - ^2., |T) , (28) 



and the second possibihty is 

(zj; S\^Si- \kl;T) 



j=l,2 



^ J2 {6,i{t\B,^\k) + 6a{3\B,^\k)-6.,{j\B,^\l)-6j,{t\B,Al)} (29) 

Xi=x,y,z 

X (5| - 52., |T) . (30) 



Matrix elements for the \T) < > \T) coupling 

The only possibihty is, 

{^j{ru rs); T| |^j(r'i, rs); T) (31) 

(u;r|^5,-Sjv|fc/;r) = 

i=l,2 

^ J2 {6,i{t\B.,^\k)-5u{j\B,^\k) + 6,,{j\B,^\l)-5,,{t\B,^\l)} (32) 



2 



X {T\ S,^^ + S2X, \T) . (33) 

The spin-couphng elements (51 Si^^ — S2xi \T) and (T| Si^^ + S2x, \T) for Xi G {x, y, z} are 
calculated straightforward using the properties of the spin operators. Their values are listed 



below for the different cases numbered from a to j. The spatial matrix elements, Eq. (27) 



(29) and (32), are composed of simple, separable integrals over the Hermite basisfunctions. 
Since 5 = over the left dot the integration in the x-direction runs over half the interval, 

poo poo 

/ da; / dy H^.A^)H^y{y)HmMH.nM^'^'^''^'''^- (34) 

Jo J-oo 
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As the Hermite functions have well defined parity we can use the values of tabulated integrals 
over the whole interval. Denoting the matrices made up by the spatial integrals above, by 
S, S — T and T respectively, we can set up the following matrix. 





IS) 


IT,) IT) IT,) 


IS) 


a 


h . c d 

$-1 


|T„) 


,b* 


e \ /I g 


IT.) 


* 

c 


/[■■| 1 i 


IT,) 


f 


g iMl\ J 



a 
b 



V2 



{B^ - iBy 



1 

71 



[B^ + iBy 



9 
h 



1 

71 
1 

71 

-B, 


B, 



{B^ + iBy) 

{B^ - iBy) 



In the Hermite basis each square in the matrix represents a ~ 4000 x 4000 matrix. Again 
we convert to the adiabatic eigenfunction basis at each electric field strength in order to 
keep the total size of the matrix small (4n x An, n = 50). 
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FIG. 5: Probabilities for the transitions to the charge localized state during the different ap- 
proaches: (second panel) using a sequence of intuitive pulses [Sj and (last panel) optimizing the 
field using the J a functional. Populations are plotted for the ground state (black curves), the sec- 
ond excited state (green curves) and the two ionic states (red and blue curves). The corresponding 
pulses are shown above both panels. The white line in the upper panel is a close-up of the last 
pulse. The final times are 852 ps and 117 ps. The optimization was done using At = 0.05. 
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FIG. 6: Properties of the optimization routines using only energy penalty (red) and including 
structure penalty (blue) in the transition to the CLS. The upper panel shows the convergence of 
the yield as a function of iteration number. The middle panel shows the optimal pulses obtained 
in the two optimizations. The dashed gray curve in the background shows the initial starting field, 
0.01 sin^(7r t/T) cos(O.Olt). The bottom panel shows the spectrum of the two pulses. The final time 
is 67 ps and the timestep used is 0.5. 
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FIG. 7: Figure showing an optimized pulse for the transition from |0) to \CLS) using a constant 
initial field (lower panel). Population as a function of time is plotted in the upper panel for the 
ground state (black curve), the second excited state (green curve) and the two ionic states (red 
and blue curve). 
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FIG. 8: Figure showing optimized transition and corresponding pulse at the anticrossing in Fig. [T} 
The upper panel shows population in the two lowest energy eigenstates during the pulse, before 
and after the adiabatic switch and evolve, indicated by the gray area (see text for description). 
The same pulse is employed in both cases. Lower panel: the optimized pulse. 
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